HIResist: a database of HIV-1 resistance to broadly neutralizing antibodies.

MOTIVATION: Changing the course of the human immunodeficiency virus type I (HIV-1) pandemic is a high public health priority with approximately 39 million people currently living with HIV-1 (PLWH) and about 1.5 million new infections annually worldwide. Broadly neutralizing antibodies (bnAbs) typically target highly conserved sites on the HIV-1 envelope glycoproteins (Envs), which mediate viral entry, and block the infection of diverse HIV-1 strains. But different mechanisms of HIV-1 resistance to bnAbs prevent robust application of bnAbs for therapeutic and preventive interventions. RESULTS: Here we report the development of a new database that provides data and computational tools to aid the discovery of resistant features and may assist in analysis of HIV-1 resistance to bnAbs. Bioinformatic tools allow identification of specific patterns in Env sequences of resistant strains and development of strategies to elucidate the mechanisms of HIV-1 escape; comparison of resistant and sensitive HIV-1 strains for each bnAb; identification of resistance and sensitivity signatures associated with specific bnAbs or groups of bnAbs; and visualization of antibody pairs on cross-sensitivity plots. The database has been designed with a particular focus on user-friendly and interactive interface. Our database is a valuable resource for the scientific community and provides opportunities to investigate patterns of HIV-1 resistance and to develop new approaches aimed to overcome HIV-1 resistance to bnAbs. AVAILABILITY AND IMPLEMENTATION: HIResist is freely available at https://hiresist.ahc.umn.edu/.

Anti-viral immune responses in a primitive lung: characterization and expression analysis of interferon-inducible immunoproteasome subunits LMP2, LMP7 and MECL-1 in a sarcopterygian fish, the Nigerian spotted lungfish (Protopterus dolloi).

Lungfishes (Dipnoi) represent the closest ancestor of tetrapods. Dipnoi have dual breathing modes extracting oxygen from water and air. The primitive lungs of lungfishes are exposed to external antigens including viruses. To date, the immune response of lungfishes against viruses has not been investigated. During viral immune responses, cell exposure to type I interferon induces the replacement of the constitutive proteasome with LMP2, LMP7 and MECL-1 beta subunits forming the immunoproteasome and enhancing antigen presentation to MHC class I molecules. In order to study the immune defense system of the lungfish lung, we have characterized for the first time the three immunoproteasome subunits in the sarcopterygian fish, the Nigerian spotted lungfish (Protopterus dolloi). LMP2, LMP7 and MECL-1 were identified in P. dolloi and their sequences encoded predicted proteins of 216, 275 and 278 amino acids, respectively. The mRNA of these three genes was expressed in multiple tissues, including the lung, with the highest abundance observed in kidney and post-pyloric spleen. In vitro stimulation of lungfish lung and kidney primary cell cultures with PolyI:C for 4 and 12 h resulted in increased LMP2, LMP7 and MECL-1 expression in both tissues. These results suggest a central role of these genes in the activation of an antiviral immune response in lungfish. Importantly, they indicate that the primitive lung of the common ancestor of all tetrapods is capable of inducing the expression of these genes in response to viral stimulation.

Conformational analysis of piperazine and piperidine analogs of GBR 12909: stochastic approach to evaluating the effects of force fields and solvent.

Analogs of the flexible dopamine reuptake inhibitor, GBR 12909 (1), may have potential utility in the treatment of cocaine abuse. As a first step in the 3D-QSAR modeling of the dopamine transporter (DAT)/serotonin transporter (SERT) selectivity of these compounds, we carried out conformational analyses of two analogs of 1: a piperazine (2) and a related piperidine (3). Ensembles of conformers consisting of local minima on the potential energy surface of the molecule were generated in the vacuum phase and in implicit solvent by random search conformational analysis using the Tripos and MMFF94 force fields. Some differences were noted in the conformer populations due to differences in the treatment of the tertiary amine nitrogen and ether oxygen atom types by the force fields. The force fields also differed in their descriptions of internal rotation around the C(sp³)-O(sp³) bond proximal to the bisphenyl moiety. Molecular orbital calculations at the HF/6-31G(d) and B3LYP/6-31G(d) levels of C-O internal rotation in model compound (5), designed to model the effect of the proximity of the bisphenyl group on C-O internal rotation, showed a broad region of low energy between -60° to 60° with minima at both -60° and 30° and a low rotational barrier at 0°, in closer agreement with the MMFF94 results than the Tripos results. Molecular mechanics calculations on model compound (6) showed that the MMFF94 force field was much more sensitive than the Tripos force field to the effects of the bisphenyl moiety on C-O internal rotation.

Toward Quantitative Structure-Property Relationships for Charge Transfer Rates of Polycyclic Aromatic Hydrocarbons.

Quantitative structure-property relationships (QSPRs) have been developed and assessed for predicting the reorganization energy of polycyclic aromatic hydrocarbons (PAHs). Preliminary QSPR models, based on a combination of molecular signature and electronic eigenvalue difference descriptors, have been trained using more than 200 PAHs. Monte Carlo cross-validation systematically improves the performance of the models through progressive reduction of the training set and selection of best performing training subsets. The final biased QSPR model yields correlation coefficients q(2) and r(2) of 0.7 and 0.8, respectively, and an estimated error in predicting reorganization energy of ±0.014 eV.

Quantitative structure-activity relationship studies of threo-methylphenidate analogs.

Complementary two-dimensional (2D) and three-dimensional (3D) Quantitative Structure-Activity Relationship (QSAR) techniques were used to derive a preliminary model for the dopamine transporter (DAT) binding affinity of 80 racemic threo-methylphenidate (MP) analogs. A novel approach based on using the atom-level E-state indices of the 14 common scaffold atoms in a sphere exclusion protocol was used to identify a test set for 2D- and 3D-QSAR model validation. Comparative Molecular Field Analysis (CoMFA) contour maps based on the structure-activity data of the training set indicate that the 2' position of the phenyl ring cannot tolerate much steric bulk and that addition of electron-withdrawing groups to the 3' or 4' positions of the phenyl ring leads to improved DAT binding affinity. In particular, the optimal substituents were found to be those whose bulk is mainly in the plane of the phenyl ring. Substituents with significant bulk above or below the plane of the ring led to decreased binding affinity. Suggested alterations to be explored in the design of new compounds are the placement at the 3' and 4' position of the phenyl ring of electron-withdrawing groups that lie chiefly in the plane of the ring, for example, halogen substituents on the 3',4'-benzo analog, 79. A complementary 2D-QSAR approach-partial least squares analysis using a reduced set of Molconn-Z descriptors-supports the CoMFA structure-activity interpretation that phenyl ring substitution is a major determinant of DAT binding affinity. The potential usefulness of the CoMFA models was demonstrated by the prediction of the binding affinity of methyl 2-(naphthalen-1-yl)-2-(piperidin-2-yl)acetate, an analog not in the original data set, to be in good agreement with the experimental value.

Novel inhibitors of anthrax edema factor.

Several pathogenic bacteria produce adenylyl cyclase toxins, such as the edema factor (EF) of Bacillus anthracis. These disturb cellular metabolism by catalyzing production of excessive amounts of the regulatory molecule cAMP. Here, a structure-based method, where a 3D-pharmacophore that fit the active site of EF was constructed from fragments, was used to identify non-nucleotide inhibitors of EF. A library of small molecule fragments was docked to the EF-active site in existing crystal structures, and those with the highest HINT scores were assembled into a 3D-pharmacophore. About 10,000 compounds, from over 2.7 million compounds in the ZINC database, had a similar molecular framework. These were ranked according to their docking scores, using methodology that was shown to achieve maximum accuracy (i.e., how well the docked position matched the experimentally determined site for ATP analogues in crystal structures of the complex). Finally, 19 diverse compounds with the best AutoDock binding/docking scores were assayed in a cell-based assay for their ability to reduce cAMP secretion induced by EF. Four of the test compounds, from different structural groups, inhibited in the low micromolar range. One of these has a core structure common to phosphatase inhibitors previously identified by high-throughput assays of a diversity library. Thus, the fragment-based pharmacophore identified a small number of diverse compounds for assay, and greatly enhanced the selection process of advanced lead compounds for combinatorial design.

Genome scale enzyme-metabolite and drug-target interaction predictions using the signature molecular descriptor.

MOTIVATION: Identifying protein enzymatic or pharmacological activities are important areas of research in biology and chemistry. Biological and chemical databases are increasingly being populated with linkages between protein sequences and chemical structures. There is now sufficient information to apply machine-learning techniques to predict interactions between chemicals and proteins at a genome scale. Current machine-learning techniques use as input either protein sequences and structures or chemical information. We propose here a method to infer protein-chemical interactions using heterogeneous input consisting of both protein sequence and chemical information. RESULTS: Our method relies on expressing proteins and chemicals with a common cheminformatics representation. We demonstrate our approach by predicting whether proteins can catalyze reactions not present in training sets. We also predict whether a given drug can bind a target, in the absence of prior binding information for that drug and target. Such predictions cannot be made with current machine-learning techniques requiring binding information for individual reactions or individual targets.

Understanding virulence mechanisms in M. tuberculosis infection via a circuit-based simulation framework.

Tuberculosis (TB), caused by the bacterium Mycobacterium tuberculosis (Mtb), is a growing international health crisis. Mtb is able to persist in host tissues in a non-replicating persistent (NRP) or latent state. This presents a challenge in the treatment of TB. Latent TB can re-activate in 10% of individuals with normal immune systems, higher for those with compromised immune systems. A quantitative understanding of latency-associated virulence mechanisms may help researchers develop more effective methods to battle the spread and reduce TB associated fatalities. Leveraging BioXyce's ability to simulate whole-cell and multi-cellular systems we are developing a circuit-based framework to investigate the impact of pathogenicity-associated pathways on the latency/reactivation phase of tuberculosis infection. We discuss efforts to simulate metabolic pathways that potentially impact the ability of Mtb to persist within host immune cells. We demonstrate how simulation studies can provide insight regarding the efficacy of potential anti-TB agents on biological networks critical to Mtb pathogenicity using a systems chemical biology approach.

Feature extraction using molecular planes for fuzzy relational clustering of a flexible dopamine reuptake inhibitor.

Six rigid-body parameters (Shift, Slide, Rise, Tilt, Roll, Twist) are commonly used to describe the relative displacement and orientation of successive base pairs in a nucleic acid structure. The present work adapts this approach to describe the relative displacement and orientation of any two planes in an arbitrary molecule-specifically, planes which contain important pharmacophore elements. Relevant code from the 3DNA software package (Nucleic Acids Res. 2003, 31, 5108-5121) was generalized to treat molecular fragments other than DNA bases as input for the calculation of the corresponding rigid-body (or "planes") parameters. These parameters were used to construct feature vectors for a fuzzy relational clustering study of over 700 conformations of a flexible analogue of the dopamine reuptake inhibitor, GBR 12909. Several cluster validity measures were used to determine the optimal number of clusters. Translational (Shift, Slide, Rise) rather than rotational (Tilt, Roll, Twist) features dominate clustering based on planes that are relatively far apart, whereas both types of features are important to clustering when the pair of planes are close by. This approach was able to classify the data set of molecular conformations into groups and to identify representative conformers for use as template conformers in future Comparative Molecular Field Analysis studies of GBR 12909 analogues. The advantage of using the planes parameters, rather than the combination of atomic coordinates and angles between molecular planes used in our previous fuzzy relational clustering of the same data set (J. Chem. Inf. Model. 2005, 45, 610-623), is that the present clustering results are independent of molecular superposition and the technique is able to identify clusters in the molecule considered as a whole. This approach is easily generalizable to any two planes in any molecule.

Aplysia temptin - the 'glue' in the water-borne attractin pheromone complex.

Temptin, a component of the complex of water-borne protein pheromones that stimulate attraction and mating behavior in the marine mollusk Aplysia, has sequence homology to the epidermal growth factor (EGF)-like domains of higher organisms that mediate protein-cell surface contact during fertilization and blood coagulation. In this work, recombinant temptin for structural and functional studies was produced in Escherichia coli using a cold shock promoter and purified by RP-HPLC. CD spectra confirmed a predominantly beta-sheet structure. Two disulfide bonds were determined via limited proteolysis and MS. One internal disulfide (Cys57-Cys77) was predicted from initial alignments with class I EGF-like domains; the second, between Cys18 and Cys103, could protect temptin against proteolysis in seawater and stabilize its interacting surface. A three-dimensional model of temptin was prepared with our MPACK suite, based on the Ca(2+)-binding, EGF-like domain of the extracellular matrix protein fibrillin. Two temptin residues, Trp52 and Trp79, which align with cysteine residues conserved in fibrillins, lie adjacent to and could stabilize the disulfide bonds and a proposed metal-binding loop. The water-borne pheromone attractin in egg cordon eluates is complexed with other proteins. Docking results with our model and the NMR structure of attractin suggest that one face of temptin interacts with the pheromone, perhaps controlling its access to the cellular receptors. Gel shifts confirmed that temptin complexes with wild-type attractin. These results indicate that temptin, analogous to the role of fibrillin in controlling transforming growth factor-beta concentration, modulates pheromone signaling by direct binding to attractin.

Flavitrack: an annotated database of flavivirus sequences.

MOTIVATION: Properly annotated sequence data for flaviviruses, which cause diseases, such as tick-borne encephalitis (TBE), dengue fever (DF), West Nile (WN) and yellow fever (YF), can aid in the design of antiviral drugs and vaccines to prevent their spread. Flavitrack was designed to help identify conserved sequence motifs, interpret mutational and structural data and track evolution of phenotypic properties. SUMMARY: Flavitrack contains over 590 complete flavivirus genome/protein sequences and information on known mutations and literature references. Each sequence has been manually annotated according to its date and place of isolation, phenotype and lethality. Internal tools are provided to rapidly determine relationships between viruses in Flavitrack and sequences provided by the user.

Recombinant production and structural studies of the Aplysia water-borne protein pheromone enticin indicates it has a novel disulfide stabilized fold.

Enticin is one of three Aplysia proteins released during egg laying that act in concert with the pheromone attractin to attract other Aplysia and stimulate mating behavior. Whereas the enticin cDNA predicts a 69-residue mature protein, enticin isolated from the albumen gland was found to be posttranslationally processed in vivo by cleavage at Arg(50) residue to generate a smaller 49-residue mature peptide. The Arg(50) cleavage site is conserved in enticin from both Aplysia californica and Aplysia brasiliana. In order to generate sufficient enticin for structural studies, recombinant full-length protein was produced in a soluble form in Escherichia coli using a cold shock promoter-based protein expression system. The enticin cDNA was cloned into the bacterial vector pCold III, and efficiently expressed, as determined by amino acid microsequence and immunoblot analyses. Recombinant enticin, which contained an additional N-terminal 13-residue translation-enhancing element, was purified by reversed-phase HPLC and compared to enticin isolated from the albumen gland. The three disulfide bonds in enticin were characterized by endoproteinase Glu-C proteolysis followed by mass spectrometric characterization of the fragments. The cysteine pairing, for both recombinant and native enticin, was I-II, III-IV, and V-VI, confirming that the protein produced in the bacterial system was correctly folded. The circular dichroism spectrum of the recombinant protein indicated it was predominantly alpha-helical. While this was consistent with fold recognition server results indicating a fold for enticin similar to that of attractin, the disulfide bonding pattern differs. A model for enticin was prepared based on its helical structure and these disulfide constraints.

Singular value decomposition of torsional angles of analogs of the dopamine reuptake inhibitor GBR 12909.

Analysis of large, flexible molecules, such as the dopamine reuptake inhibitor GBR 12909 (1), is complicated by the fact that they can take on a wide range of closely related conformations. The first step in the analysis is to classify the conformers into groups. Here, Singular Value Decomposition (SVD) was used to group conformations of GBR 12909 analogs by the similarity of their nonring torsional angles. The significance of the present work, the first application of SVD to the analysis of very flexible molecules, lies in the development of a novel scaling technique for circular data and in the grouping of molecular conformations using a technique that is independent of molecular alignment. Over 700 conformers each of a piperazine (2) and piperidine (3) analog of 1 were studied. Analysis of the score and loading plots showed that the conformers of 2 separate into three large groups due to torsional angles on the naphthalene side of the molecule, whereas those of 3 separate into nine groups due to torsional angles on the bisphenyl side of the molecule. These differences are due to nitrogen inversion at the unprotonated piperazinyl nitrogen of 2, which results in a different ensemble of conformers than those of 3, where no inversion is possible at the corresponding piperidinyl carbon.

Novel feature extraction technique for fuzzy relational clustering of a flexible dopamine reuptake inhibitor.

This paper describes a novel clustering methodology for classifying over 700 conformations of a flexible analogue of GBR 12909, a dopamine reuptake inhibitor that has completed phase I clinical trials as a treatment for cocaine abuse. The major aspect of the clustering methodology includes an efficient data-conditioning scheme where a systematic feature extraction procedure based on the structural properties of the molecule was used to reduce the associated feature space. This allowed region-specific clustering that focused on individual pharmacophore elements of the molecule. For clustering of the reduced feature set, the fuzzy clustering partitional method was utilized. Due to the relational nature of the feature data, fuzzy relational clustering was employed, and it successfully detected natural groups defined by rotational minima around N(sp(3))-C(sp(3)), O(sp(3))-C(sp(3)), and C(sp(3))-C(sp(2)) bonds. The proposed clustering methodology also employed several cluster validity measures, which corroborated the partitions produced by the clustering technique and agreed with the results of hierarchical clustering using the XCluster program. Representative structures which exhibited a reasonable spread of energies and showed good spatial coverage of the conformational space were identified for use as putative bioactive conformations in a future Comparative Molecular Field Analysis of GBR 12909 analogues. The clustering methodology developed here is capable of handling other computational chemistry problems, and the feature extraction technique can be easily generalized to other molecules.

Conformational analysis of methylphenidate: comparison of molecular orbital and molecular mechanics methods.

Methylphenidate (MP) binds to the cocaine binding site on the dopamine transporter and inhibits reuptake of dopamine, but does not appear to have the same abuse potential as cocaine. This study, part of a comprehensive effort to identify a drug treatment for cocaine abuse, investigates the effect of choice of calculation technique and of solvent model on the conformational potential energy surface (PES) of MP and a rigid methylphenidate (RMP) analogue which exhibits the same dopamine transporter binding affinity as MP. Conformational analysis was carried out by the AM1 and AM1/SM5.4 semiempirical molecular orbital methods, a molecular mechanics method (Tripos force field with the dielectric set equal to that of vacuum or water) and the HF/6-31G* molecular orbital method in vacuum phase. Although all three methods differ somewhat in the local details of the PES, the general trends are the same for neutral and protonated MP. In vacuum phase, protonation has a distinctive effect in decreasing the regions of space available to the local conformational minima. Solvent has little effect on the PES of the neutral molecule and tends to stabilize the protonated species. The random search (RS) conformational analysis technique using the Tripos force field was found to be capable of locating the minima found by the molecular orbital methods using systematic grid search. This suggests that the RS/Tripos force field/vacuum phase protocol is a reasonable choice for locating the local minima of MP. However, the Tripos force field gave significantly larger phenyl ring rotational barriers than the molecular orbital methods for MP and RMP. For both the neutral and protonated cases, all three methods found the phenyl ring rotational barriers for the RMP conformers/invertamers (denoted as cte, tte, and cta) to be: cte, tte > MP > cta. Solvation has negligible effect on the phenyl ring rotational barrier of RMP. The B3LYP/6-31G* density functional method was used to calculate the phenyl ring rotational barrier for neutral MP and gave results very similar to those of the HF/6-31G* method.